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SUMMARY 


The disturbances produced by a moving spherical diffusion flame are inves- 
tigated within the framework of a linearized theory. After the flame position 
and species concentrations are determined, the problem of determining the asso- 
ciated density, pressure, temperature, and velocity fields is delineated. Ex- 
plicit results for certain limiting situations are discussed. 


INTRODUCTION 


This paper deals with the acoustic disturbances produced by an unsteady 
spherical diffusion flame. The basic problem is envisaged as follows. Inside 
an initial sphere of radius ro, a mixture of oxidant (or fuel) and product 
gases exists, and outside the sphere a mixture of fuel (or oxidant) and product 
gases exists, as depicted in Figure 1. At an initial instant the spherical sur- 
face that separates the two initial mixtures disappears, and the subsequent com- 
bustion and acoustic disturbances are to be determined. We assume that the com- 
bustion of the initially unmixed oxidant and fuel species is confined to a spher 
ical Burke- Schumann flame surface. 

This problem is associated with the theory of particle and liquid-droplet 
combustion at high pressures, such as treated by Spalding (Ref. 1) and Rosner 
(Ref. 2). An applicable situation occurs when bubbles of oxidant (or fuel) are 
injected into a medium of fuel (or oxidant). The bubbles burst, or are ignited, 
and the subsequent disturbance field is to be determined. This investigation 
also pertains to the general study of unsteady diffusion flames and to the broad 
subject of spherical explosions. 

Most investigations of diffusion flames have dealt with steady flows, start 
ing with the original work of Burke and Schumann (Ref. 3) . A number of such 
problems are described in texts and review articles (Refs. 4 and 5). In par- 
ticular, detailed mathematical expositions for steady linearized flows have been 
developed by Clarke (Refs. 6 and 7); the present work is something akin to this 
framework of analysis. Unsteady diffusion flames have not been as thoroughly 
treated, Clarke and Stegen (Ref. 8) considered unsteady perturbations on a 
two-dimensional flame sheet. Rasmussen (Refs. 9 and 10) developed a linearized 
theory for the one«^dimensional motion induced by a diffusion flame. In this 
paper the related problem of a spherical diffusion flame is treated. 
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Whereas the general form of the problem is nonlinear, the limiting acoustic 
approximation is of some interest. The explicit analytic results obtained with- 
in the linearized theory contain the embryonic behavior of the more complicated 
nonlinear problem. The linearized results, besides having their own intrinsic 
value, can be used to estimate the magnitude of convective terms and deviations 
from constant-property approximations. The gas dynamic interactions of the 
pressure, temperature, and velocity fields can be examined coherently within a 
linearized framework. 


FORMULATION OF THE PROBLEM 
Preliminary Remarks 

We consider the basic combustion process that involves a mixture of oxi- 
dant (X), fuel (F), and product (P) species governed by the reaction equation 

^f 

XX+fFi^gP (1) 

t 

r 

where X j g are the stoichiometric coefficients. The corresponding mole- 

cular weights for the species are denoted by W^, Wp, and Wp. The forward and 
reverse reaction times are denoted by tf and t^-. The essence of combustion is 
that tf ^ 0 and t^. In this sense, the reaction associated with (1) moves 

only to the right, and this is the limiting situation considered here. This 
simple reaction model has been utilized previously (Refs. 5-10). In this limit 
a uniform oxidant-product mixture exists inside the initial spherical diaphragm, 
of radius ro, and a uniform fuel-product mixture outside (or vice versa). 

In the framework of a linearized theory, we begin with a basic uniform 
medium composed of entirely the product species at pressure, density, and tem- 
perature conditions denoted by p^, p^, and T^. Inside the initial sphere we 
suppose that the pressure and temperature are perturbed such that p = p^(l+Ap) 
and T = To(1+At). The corresponding oxidant and fuel mass fractions inside 
the initial sphere are denoted by = Ay snd Cp = 0. Outside the initial sphere 

the pressure and temperature are unperturbed, but the mass fractions are given 
by Cx = 0 and Cp = Ap (see Fig. 1). At a given instant, t = 0, the diaphragm 
disappears. We assume that the ambient temperature and density, Tq and Pq, 
are sufficient to initiate combustion spontaneously. Combustion of the oxidant 
and fuel occurs at an interface, and a diffusion flame is established. In the 
limit tf 0, the diffusion flame collapses to a discontinuity surface, and all 
the reactions occur entirely on this surface. On either side of the surface 
there are no reactions, but binary diffusion takes place. 

The limiting situation tf 0 produces a singular perturbation problem. 

The outer problem corresponds to the flame envisaged as the discontinuity sur- 
face. The inner problem deals with the structure of the flame and is obtained 
by matching with the outer problem. This procedure has been established pre- 
viously for steady-flow problems (Refs. 5-7) and for one -dimensional unsteady 
flow problems (Refs. 9 and 10). Here we restrict outselves to the lowest-order 
outer problem, with the flame treated as a discontinuity surface. 
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Basic Equations 


VThen the mass fractions Cx and Cp are small, the pressure and thermal diffu 
Sion coefficients are small and are represented by nonlinear contributions. 
Consequently, since binary diffusion prevails on either side of the flame sheet, 
the diffusion-flux vectors are represented by Flck's law. As is usual in these 
problems, we assume the binary-diffusion coefficients are equal: Dxp = Dpp = D. 

We also assume that the components of the mixtures behave as thermally perfect 
gases. We now introduce dimensionless pressure, density, and temperature per- 
turbations p* , p*, T* defined by 

P = P^(l+P') , p = p^(Up’) , T = T^(l+T*) (2) 

We have correspondingly for the mass fractions and velocity: 



V = V/a^ 

o 

where af is the frozen speed of 
of specific heats of the product 
the dimensionless space and time 


= Y P^/P^ (4) 

o 

sound of the ambient medium and y is the ratio 
species in the ambient medium. We represent 
variables by 


r = 


r/v 


T = a. 




(5) 


.f , ■ = “f 

o o 

4'--' 

where v = (2[L + X )/p , and pi and X are the first and second coefficients 

of viscosity. °Since the pr^gblem of interest is spherically symmetric we intro- 
duce a velocity potential, v = Vcp, where V = 3/Sr (the barred space variables 
are dimensional) . The characteristic Prandtl and Schmidt numbers are defined 
as 

( 6 ) 


P = p V c /k 

r OOP o 


S 3 V /D 
c o o 


where Cp^ is the constant-pressure specific heat of the product species and 
Icq is the thermal conductivity, the subscript naught denoting the ambient state. 
With the above definitions, the linearized equations for mass, momentum, species 
energy, and thermal equation of state become 


5^ + = 0 

DpT - ^ P,, I? + P^Q R 


P=v[v^cp-|£] 




p = p + T + + cVpCp 


(7a, b) 
(8a,b) 

(9a,b) 


The primes denoting the perturbation variables have been omitted. In the above 
we have the operators Dp and Dg, the parameters of Q^x and ofp, and the dimen- 
sionless heat of reaction Q defined as 


or s - 1 
X Rp 


»s ■ ^ I? - 

F Rp 


(10) 

(11) 
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Q < Cx w (h^ - hp ) + tWpdv ■*'p 

O O 0 0^0 


( 12 ) 


where Rj^, Rp and Rp are the specific gas constants. The symbol R denotes the 
dimensionless reaction rate. For the flame discontinuity sheet, R is given by 

(13) 


R = S(T) 6(Z^) 


where S(T) is the flame strength and 6(Zj^) is the Dirac delta function with Z 
the coordinate measured normal from the flame sheet. Derivations of these ^ 
equations can be found in references 6, 9 and 10, and also in references 11 
and 12 when there are no chemical reactions. 

Continuity Conditions 


All the physical flow variables are continuous at the flame sheet. Inte- 
gration of equations (8a, b) and (9a) across the flame sheet shows that the fol- 
lowing normal derivatives are discontinuous across the flame sheet: 

= S(T) (14) 

where the bracket rotation denotes the jump in value across the discontinuity. 
These jump conditions suggest that the following new variables be defined that 
do have continuous normal derivatives across the flame: 





n 

“i 

1 


(jlS] 

Un 


J_ bn 

yj_ - 

Un 

\qpJJ 


* 

C = C. 


K 

- Y ^ C 


a. P Q 

T = T + 




(15a. b) 


where 




* 


is the stoichiometric fuel-oxidant ratio. With the new variables C and T 


(16) 

* 


the reaction rate can be eliminated from the set of equations. 

can be combined and equation (9a) can be recast to read 

P Q(P -S ) Sc 
Op r^^ r c^ X 




= 0 


D T* 

P Y 


p 

r 5T + X 




Equations (8a, b) 


(17a, b) 


The linearized state equation (9b) becomes 

where 


p = p + t'' ** ^ 


P Q 

T. - ^ 

E = 


S,XWx 


-“x-“f^4 


(18a, b) 


Equation (17a) is uncoupled from the other equations and can be solved 
separately, yielding Cx and Cp as well as the flame position. A single equation 
for the velocity potential can then be found by eliminating T* in equation (17b) 
by means of equation (18) and then utilizing equations (7a, b) to eliminate p 
and p . We obtain 


L ( cp ) = ( S^-Pp [ q ^. 


* ''r'’ 

F S XW, 


TT 


+ E Dp C 


where 


L(cp) s Y V cp„ 


[v^cp 


(Y+P^) 9^^ 


J r 


„2 *1 
-V cp 


(19) 


( 20 ) 
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MASS FRACTIONS AND FLAME POSITION 


Equation (17a) for C* can be solved separately. The initial conditions 
are C*(r > rQ,0) = Ap and C*(r < ro) = -IfAp. Further, C* is finite at the 
origin and infinity, and both C* and are continuous at r=rQ. The result 
is (Ref. 2): 



^ erfc 

-d+r^)2/4T^, 


'}] 


( 21 ) 


where r^^ = r/r^ and h T/(r^S^). (22), (23) 

'k 

The flame position r = rg(T) is determined from the condition C (r^,T) = 0 
but rg(T) cannot be obtained explicitly. For small T*, however, we obtain 




.) = 1 + 2 AT^ - 2t,, + 2AT 


3/2 ,, 8 ,2^ 2 5/2^ 

* - (4 + 3 A )T^^ + 0(T^' ) 


(24) 


where 


erf (A) = (T-1)/(Y+1). 


(25) 


The position of the flame with time is a function of the initial conditions 
through the parameter defined by (16). When ^=1, the initial conditions are 
stoichiometric, and the flame travels inward toward the origin, linearly with 
T when T is small. When i > 1 the initial conditions are oxidant-rich and the 
flame initially moves outward from the origin and then subsequently toward the 
origin. The case of fuel-rich initial conditions corresponds to ^ < 1, and the 
flame moves only toward the origin, the position varying like when T is small. 
This behavior is contrary to the one-dimensional (Refs, 9 and 10) problem in 
which the flame remains stationary under stoichiometric initial conditions and 
moves either to the right or left depending on whether ^ is greater than or less 
than unity. Figure 2 shows the position of the flame as a function of T* for 
^ = 1/2, 1, and 2. When ^=1, for instance, the flame reaches the origin and 
becomes extinguished when — 0.21. 

The spherical diffusion flame travels ultimately inward. Because the 
species outside the sphere is of infinite extent and the species inside the 
sphere is of finite extent, the species inside the sphere tends to be consumed. 
Consequently, the flame e/entually will move inward in order to add. a relative 
diffusion rate for the vanishing inside species and thus maintain stoichiometric 
combustion at the flame. When the flame reaches the origin, the species ori- 
ginally inside the sphere has been completely consumed and the flame becomes 
extinguished. 

The mass fractions are determined from the function C*. For r < Tg, we 
have = -C*/(YAp) and Cp=0. For r > rg, we have Cy^=0 and Cp=C*. These 

profiles are plotted in Figure 3 for Y=1 and T* < 0.21. When the flame becomes 
extinguished at the origin, the species inside the original sphere has been 
completely consumed, in this case the oxidant. For times greater than the ex- 
tinguish-time, the species outside the original sphere proceeds to establish a 
uniform state. The return to a uniform state is shown in Figure 4. 
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SOLUTION FOR P = S and E = 0 
r c 


Laplace-Transform Analysis 


With C and known, equation (19) constitutes a fifth order non-homo- 
geneous equation for cp. Because the position of the flame varies with time, 
the problem is difficult to solve. A great simplification occurs when we set 
Pr=S^ and E=0. For most gases setting ^=3^ is a good approximation. The 
approximation E=0 is valid when ratios of specific heats of all the species are 
the same. These approximations correspond to the Shvab-Zeldovich approximation 
(Ref. 4). With these approximations, the problem is very similar to the binary 
diffusion problem of Reference 12, without pressure diffusion, and we proceed 
in similar fashion. 


If cp(r,s) denotes the Laplace transform with respect to time of cp(r,T), 
the solution of the transform of L(9) = 0 is 

rA 


r cp(r < 1 + A^ sinh sinh X^r 

° Ys 

r 9(r > r^,s) = exp(-Xj^r) + exp(-X^r) 

where A-,, A^j B- and are arbitrary constants of integration, and 
i 2 i 2 I 2 • 

^ r(Y+Pr>®^ + Pj-^s J[(Y-Pj.)s + pj + 4(P^-l)P^s^^ 

^1,2 ~ L 2(1+Ys) J 

where the plus sign holds for Xj^ and the negative sign for X 2 . The Laplace 
transforms of the pressure and density are found from (7a, b): 


(26a) 

(26b) 


(27) 


P(r,s) = Y C (i^)j.r ' ® 

r cp]/r 

(28) 

p(r < r^,s) = -(Ap/s) - 


(29a) 

pT(r > r^,s) = -(o;pAp)/s 

- (r?)j.j./(rs) 

(29b) 


The reduced temperature, T , is determined by means of equation (18). The 
velocity is determined by v = 

The constants A^, A 29 and B 2 arc determined by requiring that P, p, v, 

and be continuous at r=r . We find that 

-X r -X^r 




1 o 




*2 ■ 


2 o 




B, = (sinh X,r - X,r cosh X,r )A,/X, 
1 lo lo lo'll 

B 2 = ^^2^ ■ ^2’^o 


where 


Al = [(A + apAp)(xLs)/s + A x2/(ys2)]/(xLx^) 


(30a, b) 
(31a) 
(31b) 
(32a) 
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(32b) 

(33) 




A = ^ + /i + O' A 

p P T XX 

The inversion of the transformed variables now remains* 

Constant-Pressure Solution 

When = Sj, = 1 we have the simplified results = /s and X 2 =s// 1-pys , 
Further when Ap = 0, the factor A 2 vanishes (hence A2-B2=0) » and the pressure 
perturbation is identically zero, p(r,T) = 0. The remaining solution is en- 
tirely diffusive: 


P F F 

v(r,T) = ^ ^ 

2r VTt7 


r -(r-r )^/4T -(r+r )^/4 t-| 

e ° + (rr^+2T)e ° ‘ 


[(rr^-2T)< 


J 


P(r,T) = Ap - (Ap + a^Ap) + (C*/Ap)]/(1 + ^) 
T*(r,T) = A^^^Li - (C*/Ap)]/(1+^) 


(34) 

(35) 

(36) 


where C*/Ap is given by equation (21), and A,j,^ = A^ + (QA^/XW^^) . Wave behav- 
ior does not occur in this special case. 


The temperature at the flame sheet is determined from expression (36) by 
setting C* = Cx = 0. We obtain 

ATj 3 T*(r^,T) = [A^ + (QA^/xWjj)]/(l + +) (37) 

The increment ATf is the adiabatic flame temperature. It does not vary with 
time. The temperature is found from equation (15b) to be 

T(r ^ r^,T) = At^ + (AT^ - A^) C*/('lfAp) (38a) 

T(r ^ r^,T) = AT^ [1 - (C*/Ap)] (38b) 

This temperature distribution is shown in Figure 5 for \|r=l and A^ = 0, 

Large-Time Behavior 

An asymptotic approximation can be obtained for large times, and the re- 
sults for density, pressure, velocity and T* are similar to those obtained in 
References 11 and 12. The major difference from these results occurs in the 
temperature distribution (as contrasted with the reduced temperature T*) . _When 
Ap^O, a gasdynamic expansion wave travels inward from the initial sphere, r^, 
cooling the gas. This wave reflects from the origin and travels outward. The 
result after this reflected wave passes the position of the initial sphere is a 
reduction in temperature by an amount (y-l)Ap/Y. This temperature residual is 
subsequently eliminated by thermal diffusion (Ref. 11). Thus for times long 
after the reflected wave passes the initial sphere, the net effect is to’ change 
the initial temperature increment such that 

A^- A^ - (V-l)Ap/v (39) 

For large Reynolds numbers (Tq 1) , the wave processes are much faster than 
the diffusion processes, and under these circumstances formulas (36-38) hold 
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for the temperature long after the reflected wave passes the initial sphere, 
with the modification (39) noted. When ^ 1, the time is also modified such 
that T T/Pj. in the function C*. 

When Ap = 0, but Pj. 1, there is a weak pressure wave generated that is 
proportial to (Pr-1) (Refs. 10 and 13). This wave does not affect the temper- 
ature distribution near the flame since A<ji is left unaltered. 

CONCLUDING REMARKS 

The present analysis demonstrates the main features associated with un- 
steady spherical diffusion flames, it would be useful to continue the analysis 

without the approximations Pj- = Sc and E = 0. For the one-dimensional problem 

(Ref. 10) there are weak waves generated by the combustion at a stationary 
diffusion flame that are proportional to E. It would be interesting to inves- 
tigate these waves for a moving spherical diffusion flame. The inner problem 
associated with the structure of the flame is also worthy of further attention. 
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Figure 5.- Temperature distribution for constant-pressure solution. 
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